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Abstract. We study Mixed Mode Oscillations (MMOs) in systems of two weakly coupled 
slow/fast oscillators. We focus on the existence and properties of a folded singularity 
called FSN 11 that allows the emergence of MMOs in the presence of a suitable global 
return mechanism. As FSN 11 corresponds to a transcritical bifurcation for a desingularized 
reduced system, we prove that, under certain non-degeneracy conditions, such a transcritical 
bifurcation exists. We then apply this result to the case of two coupled systems of FitzHugh- 
Nagumo type. This leads to a non trivial condition on the coupling that enables the existence 
of MMOs. 
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1. Introduction 

Canard phenomenon, known also as canard explosion is a transition from a small, Hopf type 
oscillation to a relaxation oscillation, occurring upon variation of a parameter. This transition 
has been first found in the context of the van der Pol equation [_], and soon after in numerous 
models of phenomena occurring in engineering and in chemical reactions [4]. A common 
feature of all these models is the presence of time scale separation (one slow, one fast variable). 
A particular feature of canard explosion is that it takes place in a very small parameter interval. 
For the van der Pol system, where the ratio of the timescales is given by a small parameter 
e, the width of this parameter interval can be estimated by 0(exp(— c/e:), where c > is 
a fixed positive constant. The transition occurs through a sequence of canard cycles whose 
striking feature is that they contain segments following closely a stable slow manifold and 
subsequently an unstable slow manifold. 

The work on canard explosion led to investigations of slow/fast systems in three 
dimensions, with two slow and one fast variables. In the context of these systems the concept 
of a canard solution or simply canard has been introduced, as a solution passing from a stable 
to an unstable slow manifold [1, 12, 16]. Canards arise near so called folded singularities 
of which the best known is folded node [1, 16] and [18]. Unlike in systems with one 
slow variable, canards occur robustly (in 0(1) parameter regions) in systems with two slow 
variables. Related to canards are Mixed Mode Oscillations (MMOs), which are trajectories 
that combine small oscillations and large oscillations of relaxation type, both recurring in an 
altemating manner. Recently there has been a lot of interest in MMOs that arise due to the 
presence of canards, starting with the work of Milik, Szmolyan, Loeffelmann and Groeller 
[12]. The small oscillations arise during the passage of the trajectories near a folded node or a 
related folded singularity. The dynamics near the folded singularity is transient, yet recurrent: 
the trajectories return to the neighborhood of the folded singularity by way of a global return 
mechanism [^]. 

The setting of folded node combined with a global return mechanism, elucidated in 
[3], led to the explanations of MMO dynamics found in applications [15, 14, 8, 17]. A 
shortcoming of the folded node setting is the lack of connection to a Hopf bifurcation, which 
seems to play a prominent role in many MMOs. This led to the interest in another, more 
degenerate, folded singularity, known as Folded Saddle Node of type II (FSNII), originally 
introduced in [12] and recently analyzed in some detail by Krupa and Wechselberger [11]. 
Guckenheimer [10] studied a very similar problem in the parameter regime yet closer to the 
Hopf bifurcation, calling it singular Hopf bifurcation. The transition between the two settings 
was studied in [' ]. Another interesting singularity, that was mentioned in [ 11 ] and can lead to 
rich families of MMOs, without the presence of a Hopf bifurcation, is Folded Saddle Node of 
type I (FSN I). For a more comprehensive overview we refer the reader to the recent review 
article [ • ]. The theory of two slow and one fast variables has recently been generalized by 
Wechselberger [20] to arbitrary finite dimensions. 

In this article we study systems of two weakly coupled slow/fast oscillators. We assume 
that in the absence of the coupling one of the oscillators is undergoing a canard explosion 
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and the other is in general position. We show that turning on the weak coupling leads to 
the occurrence of MMOs. We focus on the very interesting case of FSN II, which, in the 
uncoupled case, corresponds to one of the oscillators undergoing a canard explosion, while 
the other is at a stable equilibrium. As elaborated in [ ] canard induced MMOs arise through 
a combined presence of a folded singularity and suitable retum mechanism. In this paper we 
focus on the existence and properties of a folded singularity leaving the study of the retum 
mechanism for future investigations. 

The article is organized as follows. We start with a background. Section 2, in which we 
explain the very standard case of MMOs in the context of two slow and one fast variables 
and subsequently the case of two slow and two fast variables with a simple fold line. This 
material is included for completeness and presented in such a way that the coupled oscillator 
case becomes a simple corollary. In Section 3, we treat the general case of coupled oscillators 
and in Section 4, we consider the example of two coupled FitzHugh-Nagumo systems. 



2. Background 

2.1. The basic case of one fast and two slow variables 
We consider a system of the form 

ex = f{x,y), (1) 
y = 9ix,y), X eR, y = {yi,y2) eR^,g = {gi,g2). (2) 

The functions / and g may depend on e and other parameters, but we suppress this dependence 
for simplicity of notation. The associated reduced system is 

= f{x,y), (3) 
y = g{x,y), x E R, y = {yi,y2) E . (4) 

Let So denote the curve defined by = f{x,y). Non-hyperbolic points correspond to the 
points on 5*0 for which fx = (we use the notation = df/dE). At such points the equation 
= f{x,y) cannot be solved for a; as a function of y. Suppose (0, 0) is a non-hyperbolic point. 
We make a non-degeneracy assumption fxxiO) 7^ 0. In order to obtain an explicit equation 
for the slow flow, we try to solve = /(x, y) for yi or ?/2- It is convenient to first change 
the variables to simplify the process of finding such a solution. We begin with a change of 
variables of the form x — )■ x + //(y), 77(0, 0) = 0. In the new variables, with additional scaling, 
/ has the form 

fix, yi, 1/2) = /(O, y) + x^ + 0{x^\y\,x^). 

We make a non-degeneracy assumption /j^(0, 0) 7^ (0,0). Without lost of generality (WLOG), 
we can assume that /^^(O, 0) 7^ 0. We can now introduce a new coordinate yi = — /(O, y), 
and immediately drop the tilde, to simplify the notation. In the new coordinates / has the 
form 

fix,yi,y2) = -yi + + 0(x^|?/|,x^). 
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The transformations we made may be just valid locally, i.e. only in a small neighborhood of 
the non-hyperbolic point. Now 5*0 is (locally) represented as a graph: 

yi = x'{l + 0{y2,x)). (5) 
Differentiating (5) we get 

1/1 = 2a;(l + 0(2/2, x))x + 0(x2). (6) 
Substituting (6) into (4) we get 

2x(l + 0(2/2, x))x = gi{x, x^l + 0(^/2, x)), 2/2) + 0{x^) (7) 

m = 5'2(a;,a;^(l + 0(2/2, x)), 2/2). (8) 

Now it is clear that (7) has a is singular at the fold as long as gi (0, 0,2/2) 7^ 0. Points on the fold 
line for which (/i(0, 0, 2/2) = are called folded singularities. We would like to understand 
the flow near such points better and to this end we apply a singular time rescaling to (7)-(8) 
by multiplying the right-hand side (RHS) by the factor 2x(l + 0(2/2, x)) and canceling it in 
(7). This leads to the following system: 

X = (/i(x,x2(l + 0(2/2, x)), 2/2) + O(x') (9) 
2/2 = 2x(l + 0(2/2, x))g2{x, x\l + 0(2/2, x)), 2/2). (10) 

We refer to (9)-(10) as the desingularized system. Note that folded singularities correspond 
to equilibria of (9)-(10) with x = 0. Further note that the trajectories of (7)-(8) and (9)-(10) 
restricted to the half plane x > differ only by time parametrization but are the same as sets. 
The trajectories in the half-plane x < are the same as sets but have the opposite direction 
of time. The flow near folded singularities is determined by the linearization of (9)-(10) at 
folded singularities, which is given by the Jacobian matrix: 

^i,x(0, 0,2/2) ^1,^2(0,0,2/2) 
2(72(0,0,2/2) 

Folded singularities are classified according to the type of the corresponding equilibrium 
of (9)-(10). Canards arise near folded saddles and folded nodes and small oscillations are 
associated with folded nodes, see the work of Benoit [1], and Szmolyan and Wechselberger 
[16]. Folded node, which is the only generic folded singularity whose dynamics is 
accompanied by small oscillations, occurs when the corresponding equilibrium of (9)-(10) 
has two real positive eigenvalues. 

2.2. Folded singularities 

We now consider the case when {f,g) (the RHS of (l)-(2)) depend on a regular parameter 
and (72(0,0,2/2) passes through when this parameter is varied. This passage corresponds 
to change of sign of the determinant of the Jacobian (11), from negative to positive or vice- 
versa. This, in turn, means that the eigenvalues change from one positive and one negative 
(negative determinant) to two positive (positive determinant), i.e. from saddle to node. It 
follows that this transition corresponds to the onset of small oscillations. Hence the onset of 
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small oscillations is a consequence of the passage of (72(0, 0, 1/2) through 0. This transition is 
called Folded Saddle-Node of type II (FSN II). 

As finding FSN II in coupled oscillator systems is the focus of this article w review some 
of the features of this bifurcation, referring the reader to [1 1] for details. An important feature 
of FSN II is that it corresponds to the passage of a true equilibrium of (3)-(4) through the fold 
line. We assume WLOG that FSN II corresponds to the point {x, y, A) = (0, 0, 0), where A 
is the regular parameter. Suppose that FSN II is non- degenerate, that is {x, y, A) = (0, 0, 0) 
corresponds to a transcritical bifurcation of (9)-(10). Then, for A = 0, there is a folded 
singularity satisfying (72(0, 0, 0) = and, for A 7^ but close to there is a point close to the 
origin, such that gi{x, + O(x^), ^2) = g2ix, + O(x^), 2/2) = 0, i.e. a true equilibrium of 
(3)-(4) near the fold line. For the original system (l)-(2) one can prove that this corresponds 
to a so called singular Hopf bifurcation [10], and, in a different regime of the parameter A, 
to a delayed Hopf bifurcation [1 1]. These two bifurcations lead to the existence of small 
oscillations and thus enable the existence of MMOs. We note here that the interesting case 
of FSN II is when the equilibrium is stable on the stable slow manifold and a saddle on the 
unstable slow manifold. 

For folded nodes, not close to either FSN II or to Folded Saddle Node of type I (FSN I), 
defined by gi y^ = 0, small oscillations are extremely small, which means that they cannot be 
seen, except with very detailed numerics, see [ ]. Hence the most interesting cases occur near 
FSN I or FSN II. 



2.3. Canonical system in two fast and two slow dimensions 

We consider a system with two fast and two slow variables 

ex = f{x,y,e), (12) 

y =g{x,y,e), X = (xi,X2) e M^ ?/ = (?/i,?/2) e M^ (13) 

/ = (/i) f2), 9 = (gi, 92)- Recall that in the fast formulation (12)-(13) has the form 

x' = f{x,y,e), (14) 

y' = eg{x,y,e). (15) 

Recall also that the reduced problem has the form 

= /(x,i/,0), (16) 

y = 9ix,y,0). (17) 

and the layer problem has the form 

x' = f{x,y,e), (18) 

y' = 0. (19) 

We consider a point on the fold curve (WLOG we assume that this point is the origin (0, 0)). 



Let $ : 



be the map defined by 



<l>{x,y) 



fi{x,y) 
f2{x,y) 
\ det{f^){x,y) J 



(20) 
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The fold curve is defined by the condition $ = 0. We consider a point on the fold line (WLOG 
we assume that this point is the origin (0, 0)). We now state a few conditions which assure that 
the fold curve is simple; naturally the first condition is that the linearization of layer system 
has a simple eigenvalue 0. More specifically, our first assumption is as follows: 

fxiO, 0) has one simple eigenvalue and one simple eigenvalue in the left half-plane. (21) 

Recall that in Section 2.1 we made additional assumptions, namely /xz(0,0) ^ and 
/y(0,0) 7^ (0,0). Here we generalize these conditions in the following way. Let {vi,V2) 
denote the nuUvector (the eigenvector of 0) and {wi,W2) the eigenvector of the negative 
eigenvalue. Our additional non-degeneracy conditions are as follows: 

either (0,0) or ^2 (0,0) is independent of ( ] (22) 

\W2 J 

K={det{fx)),,{0,0)-{vi,V2)^0. (23) 
We now have the following proposition. 

Proposition 1 We assume that (0, 0) is on the fold curve, i. e. $(0, 0) = 0, and that conditions 
(21)-(23) are satisfied. WLOG we assume that /j^^(0,0) is independent of {wi,W2)- Then 
there exists a (locally) unique solution of the equation $ = 0. Moreover, this solution can 
be expressed as {'Ipi{y2),4'2{y2),'ip3{y2),y2), where ip = {ipiyi^2,'ip3) : M — t- M'^ a smooth 
function, with ip{0) = (0, 0, 0). 

Proof To apply the implicit function theorem we verify that the 3x3 matrix M = 
$(a;i,x2,?/i)(0, 0) is non singular. We do this by finding three independent vectors in the 
image of M. First note that M{vi,V2,0) = (0,0, k) and that the first two coordinates 
of M{wi,W2,0) are a multiple of (101,102). Note also that the first two coordinates of 
M(0, 0, 1) are equal to fy-^ (0, 0). By assumption fy-^ (0, 0) is independent of (wi, W2), hence 
the vectors M(vi,V2, 0), M(wi,W2, 0) and M(0, 0, 1) are linearly independent. It follows that 
the equation $ = can be solved by the implicit function theorem, giving (xi, X2,yi) as a 
function of 1/2 • 

In the remainder of this section we will assume that the fold curve has been transformed 
to the coordinate line (0,0, 0,7/2) and that the fast system has been diagonalized. More 
specifically, we assume that / satisfies the following conditions 

/(0,0,0,y2,0) = 0, 1^(0,0,0,^2,0) = 0, 1^(0,0,0,1/2,0) = A(y2)<0, (24) 

0x1 0x2 

1^(0, 0, 0, 7/2, 0) = 1^(0, 0, 0, 7/2, 0) = 0. (25) 

Note that, since (wi, V2) = (1, 0) and (wi, W2) = (0, 1) the non degeneracy conditions (22)- 
(23) reduce to 

1^(0, 0, 0, 7/2, 0) ^ 0, ^(0, 0, 0, 7/2, 0) ^ 0. (26) 
WLOG we assume 

1^(0, 0, 0, 7/2, 0) = -1, ^(0, 0, 0, 7/2, 0) = I. 
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We can now expand / in Taylor series: 

( -yi + xl + 0{xiX2,xl) + 0{3) + 0{e) 
^^'^'^"1^ Xiy)x, + Oiy^,\\xr) + Oie). '• ^^'^ 

The defining conditions of the slow manifold are f{x, y, 0) = 0, or 

0= -2/i + x? + 0(xiX2,x2) + 0(3), (28) 

= X{y)x2 + Oiy,,\\xf) (29) 

From (29) we get, by the implicit function theorem, X2 = p{xi, y), with p{xi,y) = 0{yi, xf). 
Plugging into (28) we get the usual condition yi = xf + 0(xf). Following the approach of 
Section 2.1 we substitute yi = xl + 0(xf ) into (17) getting 

(2xi + 0{xl))xi = gi{xi,p{xi, xj + 0(x?), 2/2), a;? + 0(x?), ?/2, 0), (30) 
2/2 = g2ixi,p{xi,xl + 0{xl),y2),xl + O{xl),y2,0). (31) 

Finally we get the desingularized equation as follows: 

xi = gi{xi,p{xi,xl + 0{xl),y2),xl + O{xl),y2,0), (32) 
m = g2{xi,p{xi, xi + 0(x?), 2/2), X? + 0(x?), 2/2, 0)(2xi + 0{xl)) (33) 

Folded singularities are points (0, 1/2) with 1/2 satisfying (71(0, 0, 0, ?/2, 0) = 0, or, equivalently, 
equilibrium points of (32)-(33) on the fold line, i.e. satisfying xi = 0. The type of folded 
singularity is determined by the Jacobian of 



(0,0, 0,2/2) ^1,^2(0,0,0,2/2) 
2^72(0,0,0,2/2) 



(34) 



Folded saddle nodes (FSN I and FSN II) can arise in the context of systems with two slow 
and two fast dimensions in an analogous way as described in Section 2.2 for systems with two 
slow and one fast dimensions, and, in the same manner, correspond to the onset of MMOs. 

3. Coupled oscillator system 

3.1. Introduction of the system 

We consider a system of coupled oscillators in the following form 

£Xi = - 2/1 + Fi(xi) + uH({x, y) (35) 

ex2 = -y2 + F2{x2) + jyHi{x, y) (36) 

2^1 = Gi(xi, 2/1, ci) + z/i?i(x, 2/) (37) 

y2 = G2{x2, 2/2, C2) + lyH^ix, y) x = (xi, X2) G M^ 2/ = (2/1, ^2) G ^^(38) 

The parameters e and u are the singular parameter and the coupling parameter, respectively, 
and are considered to be small. The parameters ci,C2 controls the state of the uncoupled 
oscillators (moves the nuUclines). The coupling functions : — M^, i/"* : — 
and F : M2 ^ M2 defined by Rf = {H(, hI), by H' = {HI H^) and F = (Fi, F2), 
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respectively. We assume that yj = Fj{x) are S shaped curves. Written on the fast time scale 
(35)-(38) has the form: 

x[ = -yi + Fi(xi) + uHI{x, y) (39) 

4 = -^2 + F2{X2) + vHlix, y) (40) 

y[=eGi{xi,yi,Ci) + vH{{x,y) (41) 

y'^ = eG2{x2, 2/2, C2) + i^H^ix, y) (42) 

We now find the conditions for the existence of a simple fold curve, as in Section 2.3. Let 
$ defined as in section 2.3. First note that the critical manifold of (35)-(38) (the constraint 
manifold of (35)-(38) with e = 0), is defined by 

-y, + F,{x) + uHj{x,y) = 0, J = 1,2. (43) 

The linearization of the RHS of (35)-(36) is given by 



uRf F!,{x,) + uHi 



(44) 



K2 

'2,x-i ^ 2V-^2; -T- Vh2^x2 

where we assume that (x, y) is on the critical manifold, i.e. satisfies (43). The determinant of 
the matrix in (44) is as follows: 

F[{xr)F^{x2) + v{F[{x,)Hl,^ + F^{x2)HlJ + u'detHf. (45) 

Proposition 2 Let (xq, yo) satisfy $(xo, yo) = Ofor u = 0. We assume that 
Fl{xo^i) = and F2 (0:0,2) < as already specified in (24). 

And that, 

F{'{xo,i) ^ 0. 

Then, there exists a (locally) unique solution of the equation $ = 0. This solution can be 
expressed as (V'l (2/2), ^2(2/2), ^3(2/2), 2/2), where ip = (V^i, ■?/'2, '^/'a) : M Z5 a smooth 

function, with ■ip{y2,o) = (a^i.o, 2:2,0, 2/i,o)- Furthermore, the parametrization of the fold curve 
depends also smoothly on the parameters v. 

Proof 

We apply Proposition 1. For this we check the conditions (21), and (22)-(23). When 
-^1(^0,1) = ^'^d -^2(^0,2) 7^ 0, (21) hold (note that one could choose F{(a:o,i) 7^ and 
-^2(^0,2) = 0). We now have, for = 0, (^1,^2) = (1, 0) and (wi, W2) = (0, 1). It is now easy 
to verify that (22) holds automatically and (23) is equivalent to F{'(xo,i) 7^ 0. By continuity, 
the result remains valid in a neighborhood of z/ = 0. Note that, because the fast system does 
not depend on parameters ci, C2, the fold curve also does not. 



where, 



and. 
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3.2. Folded singularities and their nature 

Let a;]; = i>i{y2, i^), x*2 = ^2(1/2, ^) SM&yl = ^'3(^/2, We transform the fold curve to the line 
(0,0, 0,1/2), using the following change of variables: Xi = Xi — x\,X2 = X2 — x*2-,yi = Vi — yl- 
We omit the tilde to simplify the notation. System (35)- (38) becomes 

exi = -yi + fi{xi, 2/2) + jyh{{x, y) + 0{e) (46) 

ex2 = f2ix2, y2) + vh{{x, y) + 0{e) (47) 

iji = Gi{xi + xl,yi + yl,ci) + uH^{x + x* ,yi + yl,y2) 
dy* 

- /i(G2(x2 + x;, y2, C2) + uH^ix + x*,yi + y*, 7/2)) (48) 
dy2 

m =G2ix2 + x*2,y2,C2) + uH2{x + x*,yi + yl,y2), (49) 

/i(xi, 2/2) = F{{xl)x^ + F{'{xl)^xl + 0(x?), 
/2(x2,l/2) = F^{x;)x2 + F!l{xl)hl + 0{xl), 

h({x,y) = H(^^^{x* ,yl,y2)xi + H(,^^{x* ,yl,y2)x2 + H(^y^{x\yl,y2)yi + 0{\\{x,yi)\\'^) 

H{x2,y2) = Hl^^{x\yl,y2)xi + hI^^{x* ,yl,y2)x2 + Hiy^{x* ,yl,y2)yi + 0{\\{x,yi)\\' 

Our goal is to arrive at the desingularized system (32)-(33). The first step is to diagonalize the 
linear part of the fast subsystem (46)-(47). More precisely, we change coordinates so that the 
RHS of (46)-(49) is transformed to the form (27). Let 

■kL) - ('"TO 

be the eigenvectors of the Jacobian (44) at the points (0, 0, 0, 7/2) on the fold line, 
corresponding to the eigenvalues and A (1/2) < 0, respectively. The diagonalizing 
transformation has the form 

X = P{y)^^x, 

where the columns of the matrix P are the eigenvectors (50). In the new variables system 
(46)-(49) becomes (tilde is omitted) 

1 K 

exi = ( -yi + —xj + u{0{yi)) + 0{xiX2, xl) + 0(e)) (51) 

1 — vw 2 

ex2 = Hy2)x2 + iy{0{yi)) + 0(1 |x| n + 0(e) (52) 
yi = Gi{xi + X2w{y2) + x*, yi + yl, ci) + uH'{P{y2)x + x*, yi + yl, 2/2) 

- -^^\G2{x2 + Xiv{y2) +x;,?/2,C2) + i/if|(P(y2)a; + X*, ?/i + T/i , 1/2) j (53) 

y2= G2ix2 + xiv{y2)+x*2,y2,C2) + uH2{P{y2)x + x*,yi + yl,y2), (54) 
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where, 

K = F[\xl) + 0{u). (55) 

Note that the fast subsystem (51)-(52) is now as specified in (27). Applying the procedure 
described in Section 2.3 we obtain the reduced system 

{Kxi + 0{xl))xi = gi{xi, y2, Ci, Cs), (56) 
m = fl'2(a;i,i/2,C2), (57) 

where, 

giixi,y2,ci,C2) = Gi{xl + xi + w{y2)0{xl),yl + 0{xl),ci) 

+ uHl{xl + xi+ w{y2)0{xj),x; + v{y2)xi + 0{xl),yl + 0(x?), ys) 

- ^(G2{x*2 + v{y2)x, + 0{xl)) + i^H^ixl + x^ + w{y2)0{xl), 
dy2 V 

x*2 + v{y2)x, + 0{xl), yl + 0(x?), ys)) (58) 
92{xi,y2, C2) = G2ix*2 + v{y2)xi + 0{xl)) + uH^ixl +xi+ w{y2)0{xl), 
x*2 + v{y2)xi + 0{xl), yl + 0(x?), ^2), (59) 

and the desingularized system 

Xl = 5(l(xi,?/2,Ci,C2), (60) 

2/2 = ^2(xi, 2/2, C2){2Kxi + 0{xl)). (61) 
i.i. Folded singularities of type FSN II and the existence ofMMOs 

As discussed in Section 2.2 a well known mechanism of transition to MMOs is Folded Saddle- 
Node of type II (FSN II), see [11], [19] . This is a codimension one transition, corresponding 
to the passage of the system from a parameter region with an excitable equilibrium to a 
parameter region of MMO dynamics. It can be described as follows in the context of system 
(56)-(57): as the regular parameter is varied a stable equilibrium of (56)-(57) approaches the 
fold, and, for the critical value of the regular parameter, satisfies the conditions: 

5-1(0, 2/2, ci, C2) = (62) 

92iQ.y2.c2) =0. (63) 

On the other side of criticality there are MMOs as well as an equilibrium, now unstable. 

In the following, we assume that the parameter C2 is fixed. Note that for u = yl 
does not depend on y2 and thus dyl/dy2 = O^v). We assume that, for u = 0, the equation 
G^2(a^2' 2/2, C2) = has a unique solution x^ = ^2(2/2, C2) which is a stable equilibrium of 
the uncoupled (x2, 2/2) subsystem, and that 6*2,^2 7^ at this point. By the implicit function 
theorem this gives a unique solution y2 for the equation (63), for u small enough. Recall 
that, by hypothesis, in the case z/ = we assume that, each uncoupled subsystem (39)-(41) 
and (40)-(42) admits a unique stationary point, one attractive on the stable manifold and one 
repulsive near the fold. This correspond to the nuUcline configuration shown in Figure 1 . Let 
us denote this particular value I/2. It follows that this value I/2 determine particular values of 



Weakly coupled two slow- two fast systems 



11 




Figure 1. The 1/ ^ limit configurations. Figure a) The uncoupled system (39)-(41) admits 
a unique repulsive stationary point near the fold. Figure b) The uncoupled system (40)-(42) 
admits a unique attractive stationary point on the stable manifold. 



{xl,yi,X2) = {xl,yl,X2). We assume that solving equation (62) with these values gives a 
unique value of ci = ci. 

Recall from the discussion in Section 2.2 that a non-degenerate FSN 11 singularity 
corresponds to a transcritical bifurcation for system (60)-(61). We will show that such a 
transcritical bifurcation exists under some non-degeneracy assumptions on (35)-(38). Let us 
recall the general definition of a transcritical bifurcation, see [9, 5, 13] : 

Definition 1 Let 

X = c) (64) 

be a system of differential equations in M" depending on a real parameter c. Let us assume 
that {x, c) is a sationary point of (64): X{x, c) = 0. We say that this system admits a tran- 
scritical bifurcation at (x, c) if there exists an open neighborhood U of x in M", an open 
neighborhood V of c in M, and Ai, A2 : H- M, ^1, ^2 : ^ functions of class such 

that: 



(i) DxX{x, c) admits an eigenvalue equal to zero. This eigenvalue is simple and there is no 
other eigenvalue with real part equal to zero. 

(ii) For all c G V, ^1, ^2 the only stationary points of (64) in U, furthermore we have: 

^j(c) = X and ^1 7^ ^2 for c^c. 

(Hi) For all a & V, Aj(c) is an eigenvalue of D^Xi^^i, c), with Aj(c) = and A1A2 < 0/or 
c 7^ c. 

The following proposition establishes the existence of a transcritical bifurcation for (60)-(61). 



Proposition 3 Let us assume that, for system ( 60)-( 61 ): 
dgi dGi dG2 



dy2 dxi dy2 



(0,y2,ci) < 0(resp> 0), (65) 
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then for c\ in a neighborhood of C\ = C\ and v small enough, there exists two stationary 
points for system (60)-(61): (xi^e, 2/2,e) cind (0, y2joid)- Furthermore, if we assume that 

dci dci 

then, for u sufficiently small, there is a transcritical bifurcation. It follows that (0, y2), is an 
FSNII point. 

Proof 

The folded stationary point, (0, y2joid) is obtained by solving the equation gi{0, y2, ci, C2) = 
0, i.e, gi = 0, xi = 0. For this, we apply the implicit function theorem to gi{0, 7/2, ci, C2) as 
a function of y2 and Ci. For y2 = y2 and Ci = Ci, we have gi = 0. From hypothesis (65), it 
follows that ^(0, ^2, ci, C2) 7^ 0, then the implicit function theorem gives the existence of a 
stationary point (xi = 0, 2/2 = y2joid) for ci in a neighborhood of Ci. The second stationary 
point, (xi_e, y2,e), is obtained by solving the equation gi{xi,y2, Ci) = 0, (72(3^17 1/2) = 0. Note 
that this second stationary point corresponds to a true stationary point of system (35)-(38). 
For z/ = 0, we know that these equations admit a unique solution. Using hypothesis (65), 
we can apply the implicit function in a neighborhood of u = 0. This gives the existence of a 
unique stationary point for u sufficiently small. For the first stationary point, the folded one 
(0, y2,foid) the jacobian is given by 

J = ^ ^ ' dy2 \ (57) 

\ 2Kg2 J 

This gives: 

detJ = {^)2Kg2. 
dy2 

For the second stationary point, the ordinary one, the jacobian is given by: 

G,,,, + 0{u) + 0{x,) 
KG2,x,xi + 0{xl) KG2,y^xi + 0(x?) 

This gives: 

detJ = KGi^xiG2,y2Xi + O^xiv). 

From hypothesis (66), it follows that system (60)-(61) admits a transcritical bifurcation at 
point (0, ^2) when ci crosses the value Ci. 



(68) 



4. Example - coupled FitzHugh-Nagumo systems. 

4.1. Simple fold line 

We consider the following system 

exi = F{xi) - yi + ai{x2- xi) 
ex2 = F{x2) -y2 + a2{xi - X2) 
iji = xi - ci + I3i{y2 - yi) 

2/2 = X2 - C2 + /32(2/l - 2/2) 
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with 

F{z) = -z^ + 32. 

We define 

l{xi,yuX2,y2) = < ^ , . ^ (70) 

y F{X2) -y2 + Oi2{Xl - X2) 

Points on the fold line are defined by the condition 

det{Dl{x,y)) = 

h{x,y) = (71) 
l2{x,y) = 

wich, for system (69), is given by 

{F'{xi) - ai){F'{x2) - 02) - aia2 = 

yi = F{xi) + ai{x2 - xi) 
y2 = F{x2) + a2{xi - X2). 

Thanks to proposition 2, the existence of a simple fold line follows from F'{xi) = 0, 
F'{x2) ^ and F"{xi) ^ (the alternative choice would be F'{x2) = 0, F'(xi) ^ and 
F"{x2) 7^ 0). We obtain a solution of system (4.1) with (a:i,X2,2/i) given as function of 
y2. We denote this solution by {xl.xl.yD = {^i{y2),^2{y2),^z{y2))- Note that we have 
obtained a curve in M'' parametrized par y2, see 2 for a general statement. 



4.2. Folded singularities and their nature 

System (46)-(49) in context of (69) becomes 

' exi = -yi + fi{xi,y2) + ai{x2 - xi) + 0{e) 

£^2 = /2(a;i,y2) + «2(xi - X2) + 0(e) 
< yi = xl + xi - ci + (3i{y2 - {yt + yi)) (72) 

+ X2-C2 + ^2iyl + yi- y2)) 
, y2 = X2 + X2 - C2 + I32iyl + yi - y2) 

where 

f,(x^)=F'{xl)x, + F"{xl)^-xl 

and 

f2{x2) = F'{x;)x2 + F"{x;)^-xl 

As in section 3.2, we diagonalize the linear part of the fast system. The eigenvalues of 

Df{xl,x*2)wce: 

and A(|/2) = F'{xl) + F'{xl) - ai - 02- 

As explained in section (3.3), we focus on dynamics decribed by configuration in Figure 1, 
thus F'{xl) < 0. The associated eigenvectors are, respectively. 
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Pi 



1 

012 



F'{x*)-a2 



and Po 



a2-F'{x*) 
1 



We now diagonalize the fast subsystem using the transformation 



X = PX with X 



Xi 
X2 



where P is the matrix whose columns are the eigenvectors. We introduce the following 
notations: 

Oil a2 , 1 



w 



V 



d 



After some transformations , system (72) becomes: 



1 — vw 



exi = d{{F"{xl) 
-w{F"{x 



*\ {X1+WX2)'^ 



[Xi + WX2] 



- [bxi + X2f)) + 0(e) 



ex2 



yi 



y2 



My2)x2 + d{-v{F"{xlY^^±^ - (Xi + WX2f 

+F"(^*) (Hl+£i2)! _ (^a;^ + j:^f) + 0(e) 

xl + xi + WX2 - ci + /3i(?/2 - (yl + yi)) 
~£(^2 + bxi + X2-C2 + ^2{y*i +yi- z/2)) 
X2 + vxi +X2-C2 + (32{yl +yi- y2)- 



Hence, the reduced equation is given by 




rf((F"(xt)^^i±^^ - (Xi + WX2) 



yi] 



w{F'\xl)^^^^^^^±^-{vx, + X2f)) 



\{y2)x2 + d{-v{F"{xl 
+ F"{xl)^^^^^^±^- 



{Xi + WX2) 



[yxi + ^2)^) 



(xi + wx2f - yi) 



yi 



xl + xi+ WX2 - ci + /3i{y2 - {yl + yi)) 

dyl 



dy'. 



■{x*2 + X2-C2 + hivl + yi- y2)) 



y2 = x*2 + vxi + X2-C2 + (32{yl + yi- 2/2) 



(73) 



(74) 



(75) 



(76) 



Now, using equation (74) and the implicit function theorem in (75), as \{y2) < 0, one can 
obtain X2 as a function of Xi, 2/2 wich leads to, 

X2 = 0{xl). 

Now, (74) can be rewritten in the form 

y^ = ^xl + 0{xl) (77) 
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with 

K = F"{xl) - wv^F"{xl) 
= F"{xl) + 0{u^), 

where < z/ < max(ai, 0^2, P2)- Subsequently, we derivate (77), and plug into (75) and 
(76). We obtain, 

{Kxi + 0{xl))xi = X* + xi - ci + /3i(?/2 - yl) 

rlii* 

[xl + hxi - C2 + (32{yl - y2)) + 0{xl) 



dy2 

y2 = x*2 + bxi - C2 + /32iyl - y2) + 0{xl) 

This gives the following desingularized system: 

xi = xl + xi - ci + f3i{y2 - yl) (78) 

- ^{x; + hx, -C2 + Hyl - y2)) + 0{xl) (79) 
dy2 

= {x; + bxi -C2 + (32{yl - y2) + 0{xl)){Kxi + 0{xl)) (80) 

(81) 

Hence using notations of section 3.3, we have: 

dv* 

9i{xi, y2) = x\ + xi-ci + (3i{y2 - yl) - ^{x*2 + hxi - C2 + (32{yl - ^2)) + 0{xl) 

oy2 

and 

g2{.xi, y2) = xl + bxi -C2 + I32{yl - 2/2) + 0{xl). 

Proposition 4 Let us assume C2 « —1 and ai < f3i. For ai, ^2, (^2 small enough, let 
1/2 be the solution 0/(72(0, 1/2) = and define ci = + /3i(y2 — yl)- Then, for ai, 02, /32 
in a a small neighborhood of zero, the system (80)-(81), admits a transcritical bifurcation at 
point (xi = 0, ?/2 = 1/2, Ci = Ci. As the parameter ci crosses the value Ci from left to right, 
the folded stationary point passes from a folded saddle to a repulsive folded node, whereas 
the ordinary stationary point passes from a repulsive node to a saddle. 

Proof 

This proposition is a specific case of proposition (3). We will verify that: 

^(0,y2,ci) >0, (82) 



dy2 
dGi dG2 
dxi dy2 

dxi 



< 0, (83) 



> 0, (84) 
> 0. (85) 



dci 

dg2i0,y2joid{ci)) 
dci 
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We first show the existence of two distinct stationary points for system (80)-(81). We start 
with hypothesis (82). Here, this hypothesis reads as, 



It is sufficient to have: 



This holds since = C2 + 0(z/), C2 << —1 and ai < Pi. Now, we deal with hypothesis 
(83). Here, it reads: 

^ < 0. (86) 

dy2 

This holds since ^ = -pjj^ ~^ 0{u) < for sufficiently small. Hypothesis (84) is verified 
smce we have ^ + 0{v) > or sufficiently small. 

Finally, we come to hypothesis (85). From (82), we know that ^(0, y2, Ci) > 0. From 
(83), we know that J| < 0. Since we have, < 0, it follows that < 0, and then < 0. 



5. Numerical simulations 



We proceed numerical integration of system (69) on time intervall [0, 500], using a Runge- 
Kutta 4 method with time step of 0.001. The parameter values are the following e = 
0.1, ai = (3i = a2 = (^2 = 0.01, C2 = —1.5. The parameter ci is varied in a small 
neighborhood arround the value —0.99 . In figure 2, we represent the results of the simulation 
when ci = —0.99. Figure 2, a) represents the time-evolution of the variable Xi. Figure 
2,b) represents the behaviour of the system in the (xi,?/i) phase-plane. It shows, that the 
system evolves asymptotically toward a stable stationary point. In figure 3, we represent the 
results of the simulation when ci = —0.9883. Figure 3,a) represents the time-evolution of the 
variable Xi. Figure 3,b) represents the behaviour of the system in the (xi, yi) phase-plane. 
It shows that the system evolves asymptotically toward a limit cycle of small amplitude. In 
figure 3 , we represent the results of the simulation when ci = —0.988295. Figure 4,a) 
represents the time-evolution of the variable Xi. Figure 4,b) represents the time-evolution of 
the variable X2. Figure 4,c) represents the behaviuor of the system in the (xi, yi) phase space. 
Figure 4,d) represents the behaviour of the system in the (x2, 2/2) phase space. It shows the 
emergence of the MMO's. Figure 5,d) represents the evolution of the system in the (x2, 1/2) 
phase space. Figure 5 represents the reduced manifold: (xi, X2, yi = F(xi) + ai(x2 — xi)), 
and the behaviour of the system in the (xi, X2, yi) phase-space. It shows the emergence of the 
MMO's. 
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Figure 2. Simulation of system (69) for parameter ci = —0.99. Figure a) The time-evolution 
of xi. Figure b) The behaviour of the system in the {xi, j/i) phase-plane. It shows that the 
system evolves asymptotically toward its stable stationnary point. 



■'mmmfmm#mmmmm\ ■ 



(a) (b) 

Figure 3. Simulation of system (69) for parameter ci = —0.9883. Figure a) The time- 
evolution of xi. Figure b) The behaviour of the system in the {xi,yi) phase-plane. It shows 
that the system evolves asymptotically toward a limit cycle of small amplitude. 
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Figure 4. Simulation of system (69) for parameter ci = —988295. Figure a) The time- 
evolution of xi. Figure b) The time-evolution of xi. Figure c) The behaviour of the system in 
the {xi, yi) phase-plane. The behaviour of the system in the (x2, 2/2) phase space. It shows 
the emergence of MMO's. 




Figure 5. Simulation of system (69) for parameter ci = —988295. The reduced manifold: 
{xi,X2, 2/1 = ^"(2^1) + cti{x2 — xi) and the behaviour of the system in the {xi,X2,yi) phase- 
space (green curve). It shows the emergence of MMO's. 
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